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Abstract 

We present the electronic properties of massless Dirac fermions characterized by geometry and 
topology on a graphene sheet in this chapter. Topological effects can be elegantly illuminated 
by the Atiyah-Singer index theorem. It leads to a topological invariant under deformations on 
the Dirac operator and plays an essential role in formulating super symmetric quantum mechanics 
over twisted Dolbeault complex caused by the topological deformation of the lattice in a graphene 
system. Making use of the G index theorem and a high degree of symmetry, we study deformed 
energy eigenvalues in graphene. The Dirac fermion results in SU(4) symmetry as a high degree 
of symmetry in the noninteracting Hamiltonian of the monolayer graphene. Under the topological 
deformation the zero-energy states emerge naturally without the Zeeman splitting at the Fermi 
points in the graphene sheet. In the case of nonzero energy, the up-spin and down-spin states have 
the exact high symmetries of spin, forming the pseudospin singlet pairing. We describe the peculiar 
and unconventional quantum Hall effects of the n = Landau level in monolayer graphene on the 
basis of the G index theorem and the high degree of symmetry. 
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I. INTRODUCTION 



Carbon provides a fundamental material for all life and physical science. Carbon-based 
systems reveal a variety of structures with a great deal of physical properties. These physi- 
cal properties result from the dimensionality of the structures among systems with carbon 
atoms. For a long time, in material science, both experimentalists and theorists have sought 
for the existence of a true two-dimensional (2D) material with the thickness of a single atom, 
or a membrane of atomic thickness. This 2D material was theoretically first studied on a 
monolayer of graphite by Wallace Experimentally in the year 2004, a group led by A. K. 
Geim at the University of Manchester, U. K., realized such a 2D material under the name 
of graphene {2, 3]. 

Graphene is composed of carbon atoms placed at the vertices of a two dimensional hon- 
eycomb lattice. It is regarded as a large molecule of carbon atoms which become strongly 
bound together on the sites of the honeycomb lattice. For each carbon atom on the lattice, 
three of the four outer electrons get strongly bond with its neighboring atoms by a orbitals. 
The 2p z orbital of the fourth electron produces a ir bond with a neighboring carbon atom. 
The cr bonds form the covalent structure with a honeycomb geometry. The bond strength 
furnishes the flexibility and robustness for the lattice geometry. On the other hand, the 
71 bonds generate the intrinsic electronic structure of graphene. Each tc bond yields the 
half-filled electrons of p orbital to tunnel from a carbon atom to the neighboring one. Thus 
graphene should be regarded as a many body system on which electrons can get correlated 
from site to site, resulting in a rich collective behavior. The correlated behavior can be 
represented by quantum effects which can influence on graphene's electronic properties 

The electronic structure can be described by 2D massless relativistic fermions ^[j?! in 
graphene. The massless fermions enable us to study topological effects on electronic proper- 
ties of graphene. Topological effects are represented by the global properties of geometrical 
objects rather than their local ones. By varying the geometry, we can produce topologically 
different configurations such as a sphere or a torus on which the effective Dirac operators 
of massless fermions are well defined. Described by the Dirac fermions, graphene can have 
extraordinary properties of stability obtained in terms of geometry and topology of the un- 
derlined lattice. Apart from the robust structure of geometry, the topological properties 
can emerge due to long range quantum coherence in graphene. It follows that we can move 



2 



electrons coherently through the whole graphene molecule, resulting in its detection of ge- 
ometry or topology. This allows us to study a great deal of physical properties revealed by 
interplay between geometry and topology and quantum effects [8|, |9|. 

As quantum effects, there exist unconventional quantum Hall effects (QHE) which form a 
series of filling factors v = ±2, ±6, ±10, • • • as the four-fold degeneracy combined by spin and 
sublattice valley ones js, 7]. The energy dispersion shows a linear spectrum by the massless 
Dirac fermions with a Fermi velocity Vp ~ 10 6 m/s. The 4-fold degeneracy of the Landau level 
(LL) is lifted into 4 sublevels in the presence of an external high magnetic field. In the case of 
the tilted high magnetic field to the graphene plane, the spin degeneracy can be lifted at the 
firsr LL, resulting in the filling factor v = 4 QHE of monolayer graphene lfj. Furthermore 
for bilayer graphene, the LL spectrum is composed of eightfold degenerate states at the zero 
energy and fourfold ones at finite energies under the high magnetic field. This can allows 
us to observe the quantum Hall plateaus at a seris of v = ±4, ±8, ±12, • • • llMl3|. The 



charge carriers are chiral massive fermions which produce a parabolic energy band. The 
chiral fermions offer the unconventional integer QHE of the zero-LL anomaly which exhibits 
metallic behavior under the condition of low carrier densities and high magnetic fields in 
contrast to the conventional insulating phenomena js, 7, 9 14 1 . 

Topological configurations can produce a crucial effect on the quantum states of a system. 
In particular, they can provide the possible quantum ground states which a system can have. 
This remarkable result is described in terms of the index theorem initiated by Atiyah and 
Singer 15[. It gives the relationship between the analytic properties of the operator and 
the topological characteristic of the manifold upon which the operator is defined. The Dirac 
operator can be related to topological effects which is elegantly illuminated by the Atiyah- 
Singer index theorem in graphene [8J . It leads to a topological invariant under deformations 
on a Dirac operator and plays an essential role in formulating supersymmetric quantum 
mechanics (SUSY QM) on the graphene sheet 

In a theoretical sense, there has been at least the quantum mechanics (QM) of particles 
described by both fermionic and bosonic degrees of freedom. The SUSY QM may be hidden 
in the quantum mechanics of a particle acting on a group manifold which can be represented 
by a high degree of symmetry [17J. As an example, a spin precessing in a magnetic field 
can have the hidden SUSY. In particular, under the uniform magnetic field the LL for an 
electron can be expressed by the spectrum of the SUSY oscillator which is composed of 
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fermionic and bosonic ones. It is remarkable that this SUSY QM can possibly emerge in 
graphene with low carrier concentration and high mobility. The supersymmetry is built up 
over the Dolbeault complex due to the topological deformation on the lattice in a graphene 



system 



We exploit the G-index theorem and a high degree of symmetry to understand unusual 
quantum Hall effects of the n = Landau level in graphene. The Dirac fermion results 
in SU(4) symmetry as a high degree of symmetry in the noninteracting Hamiltonian of 
the monolayer graphene. The high symmetries in graphene sheets can not couple to an 
external magnetic field. In the absence of the magnetic field the index theorem can provide a 
relation between the zero-energy state of the graphene sheet and the topological deformation 
of the compact lattice. Under the topological deformation the zero-energy states emerge 
naturally without the Zeeman splitting at the Fermi points in the graphene sheet. In the 
case of nonzero energy, the up-spin and down-spin states have the exact high symmetries of 
spin, forming the pseudospin singlet pairing. We describe the peculiar and unconventional 
quantum Hall effects of the n = Landau level in monolayer graphene on the basis of the 

Lin 

index theorem and the high degree of symmetry [9|, |1_8[ . 

This chapter is written as follows. We explain basic properties of graphene in section II. 
In section III, we discuss a path integral of coherent states in brief. In section IV, supersym- 
metry is introduced in graphene. In the following section, we investigate the Atiyah-Singer 
index theorem and topological properties. The G index theorem and deformation is covered. 
In section VI, SUSY QM and higher spin symmetry are described. Next we illuminate the 
low energy spectrum and unconventional quantum Hall effects in monolayer graphene. And 
finally we come to summary and conclusion. 



II. BASIC PROPERTIES OF GRAPHENE 



Graphene is a molecule that is composed of carbon atoms placed on a two dimensional 
honeycomb lattice. The basic plaquette of the lattice has a hexagon and the atoms are 
located at the sites of the lattice. Electronic properties of graphene can be described by the 
tight binding model on which spinless electrons move from site to site along the links of the 
lattice without interaction each other. Under the tight-binding approximation graphene can 
be expressed by a simple Hamiltonian of coupled fermions on a hexagonal lattice {2I, [jij]. 
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The model Hamiltonian is a form given by 



H = £ <fe. (1) 

<»J> 

where < i,j > indicates nearest neighbors on the lattice. c\ and q are the creation 
and annihilation operators of the fermions located at site i with anticommuation relation 
{ j Cj } 5{j . 

In order to calculate the spectrum of Hamiltonian (1), we account for a periodicity of 
honeycomb lattice which leads to a Fourier transformation. The periodic structure provides 
the energy eigenvalue problem for the Hamiltonian in a unit cell. The unit cell consists of 
two neighboring carbon atoms called A and B. They can be expressed by the three vectors 
Ui, \/i = 1,2,3. Under the Fourier transformation of c(p) = ^2 i e ip ' Ut Ci, the Hamiltonian is 
rewritten in terms of 

^ff^,^)(. ° .„ H), (2) 

3 J J \ > \E i= i^ ip - Ui J \c B (p) J 

where ca{p) and cb(p) denote the Fourier transformed operators corresponding to the carbon 
atoms A and B, respectively. 



Now it is easy to take 
sion energy is given by 



;he eigenvalue of the energy for electrons of graphene. The disper- 
9. 



EV \ x^Fj, . q 2 V3Py . . 3?-,. 3p y 

E(p) = ±— y 1 + 3cos^ - + 4 cos — cos (3) 

where the lattice distance between atoms becomes normalized to the unity. From the 
dispersion relation obtained above, graphene can have two independent Fermi points, 
p = K± = ±^(1, i). We can expand it and then linearize it near the conical singu- 
larities of the Fermi points. Corresponding to the K + and K- at the half-filling case, the 
Hamiltonian is expressed by the Dirac operators 

H ± = ±^ F J2 7%, (4) 

H=x,y 

where p^ = —ihd^ is the covariant momentum and the Dirac matrices 7 M indicate the Pauli 
matrices 7^ = cr^. Hence the low energy theory of graphene is described by means of free 
fermions. 
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The Hamiltonian can be written in the matrix form 




H 



± 



D ± 











) 



(5) 



Here D± is a Dirac operator given by 



D± = ±v F (a x p x + a y p y ), 



% 



d z 



) 



(6) 



Here we have expressed D± in terms of d z = \{d x — id y ) and d z = \{d x + id y ). On the 
complex coordinates we can write 



So far, we have discussed the collective behavior of graphene electrons that can be governed 
by the Dirac equation. In particular, the velocity of the electrons is effectively 300 times 
smaller than the speed of light. 

Next let us describe a curved graphene. On the curved surface of graphene the low energy 
physics can be illuminated by the Dirac equation that is defined on the corresponding curved 
manifold. The curvature generates an gauge field of magnetic flux going through the curved 
graphene. This yields to a picture about how to interact gauge fields with Dirac fermions. 
In order to have the way associated with the gauge fields, let us take into account a good 
method to include curvature to graphene. The simplest way is that we cut a | piece of 
triangle from a graphene sheet and then gluing the opposite ends of the lattice. This process 
results in a single pentagon at the apex of the generated cone while all the other plaquette 
keep a hexagon. Due to the minimal geometrical distortion, the honeycomb lattice has a 
positive curvature. The curvature can be obtained by calculating a circular tangent vector, 
V, around the apex by § V ■ df = |. The generation of a single pentagon gives rise to a 
dramatic effect on the spinor, resulting in deformation of the lattice. If the spinor is parallel 
transported around the apex by an angle 2tt, it is forced at some point to make a jump 
from a site A to a site A while every site A takes only B neighborhoods or vice versa. This 
motion enables us to have the effect that the magnetic field gives on the wave function of 
a particle moving on a closed path. A full circulation provides accumulation of all phase 
factors to the particle wave function, and so generates the enclosed magnetic flux. This flux 
allows us to observe the Aharonov Bohm effect. Quantum mechanically, it can provide a 
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Ad z d- Z . 



(7) 
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discontinuity. Thus we have simultaneously to describe it through the wave function of the 
particle in the processes of being static or moving along the closed trajectory. This produces 
a vector potential term in the Hamiltonian that leaves the theory to be consistent. 

By similar procedures made on the curved graphene that compensates the jump in the 
components of the spinor, we should take into account a nonabelian vector potential A in 
the effective Hamiltonian. Around the apex we can take the circulation of A along a path. It 
is expressed by §^A-dr= |-T2, where T2 denotes the second Pauli matrix which couples the 
K + to the K- components of the spinor. The effective gauge theory can be emerged due to 
the geometric deformation on geometric variants of graphene as topological effects. In next 
section we describe a path integral of the coherent states for geometrical and topological 
properties. 



III. PATH INTEGRAL OF COHERENT STATES 



Let us construct a path integral of the coherent states on a group manifold. Following 
the Stone's approach [17] . we discuss the coherent states on a general group G. Let us take 
D(g), any g G G, as an irreducible representation of the group. Assume that |0 > is some 
state in the space of representation. Then one can define \g > by 

|<7>=£>(<7)|0>. (8) 

On the basis of the irreducible representation, Shur's lemma holds that 

1 



ncij- fXfl = 1 (9) 

where yfc) * s ^ ne vomme °f the group manifold. In the expression of Eq. (9) d[g] is the 
Haar measure on the group. 

Let us compute a thermodynamic partition function 

Z = Tr(e-? n ) (10) 

where (3 is the imaginary time. In the procedure of calculating the partition function, the 
trace is constrained to the representation space on which D(g) acts. Dividing the Matsubara 
time-interval into n parts, and using the Shur's lemma of Eq. (9), we can write down an 
iterated integral 

Tr(e^ H ) = const. J(d[g]d[g) ■■■)< g\e~^ n \g >< g\e~^ n ■ ■ ■ \g > . (11) 



7 



where const, is the constant value taken on the representation space. When taking into 
account short time intervals, one can express g ~ g + 5g, and 5g ~ 0(St), so that 

1 1 (12) 
« 1+ < ^ > + < 0|(-<JtH)|0 > +0(5t 2 ). 

On taking into consideration up to order 0(St), we can write down the formal path-integral 
expression 

TW e -0H) 

5 ( 13 ) 

= to / ^]exp(i < <?l<^ > - Jo ^ < g\R\g >)• 

Here e%] is regarded as the path-integral measure given by the Haar measure at each time 
step. 

The path-integral expression can be identified as a path integration over a quotient space 
of the group. In particular it is noted that a set of the \g > can be different from only a 
phase so that the integrand is not sensitive to the phase factor. Now suppose that H is the 
subgroup of G, constructed from exponentiating a maximal commuting set of generators, 
i.e., a maximal torus. Then |0 > is expressed by an eigenstate of the generators of H which 
means a state of definite weight. The \g > are represented by all phase multiples of one 
another in any one coset of G/H. And hence the coherent states can be described in terms 
of a bundle over G/H with the maximal torus as the gauge group, and the integration is 
made over the path in G/H . 

Let us express the integrand in the coherent-state path integral without any choice of 
representatives. In order to make a natural procedure on the independent choice of the 
representatives, we define the projection operators as 

P(g) = \gxg\. (14) 

They can be directly projected onto the physically distinct states since they do not have any 
phase ambiguity. In the integrand, the first term § < g\5g > can be taken to be a gauge 
invariant form by using Stokes theorem 

* < g\dg >= d< g\dg > . (15) 
Jr=dn Jn 

Making use of the identity form 

d < g\dg >=< dg\dg >= -Tr{dPPdP), (16) 



the first term yields 

I < g \dg >=- [ Ti(dPPdP). (17) 
Jr=an Jn 

The second term in the integrand can be rewritten by 

<g\K\g>=Tr(P(g)R). (18) 



Hence combination of Eqs. (17) and (18) leads us to the path integral given by [17] 

Tr(e^ H ) = y^rjj^ f d[g]exp (- jf Tr(dPPdP) - jf Tr(Ffo)H)) . (19) 

It is noted that the elements of G/H can only contribute to the integrand, so that the 
volume factor of the gauge group H gets out of the path integral. 

There are certain subgroups corresponded to the Lie algebra Q of a Lie group G. Let 
us describe some basic facts about semisimple Lie algebra in brief. Then it is known that 
the generators of the Q can have a decomposition into a maximally commuting set H. = 
{Hi}, i.e., the Cartan subalgebra, and a set of ladder operators, E a , one for each root vector 
a G R. The ladder operators are needed for complexifing the algebra to Q c which results 
from the group parameters to get the complex values. The Hi and the E a hold that 

[Hi, E a \ = ai E a . (20) 

Assume that |A > is an eigenvector of the Hi with eigenvalues \ such as 

Hi\X >= Xi\X > . (21) 

Then 

E a \X >= \X + a > . (22) 

The roots may be classified into two sets by an arbitrary hyperplane on a root space. On 
one set the root objects are positive roots indicated by a G R + as increasing the weights 
while the others are negative roots, by a G R + as decreasing the weights. And thus we can 
take the greatest weight as a state which is annihilated by all E a ,a G R+. 

This decomposition of the Lie algebra is made on some kinds of special subgroup. The 
Borel subgroups B± are constructed by exponentiating the algebras B± which are spanned 
by the E a , Hi, a G R±, Hi G H c . Any g G G can be decomposed into the gaussian factors 

g = C-K + ; C^eZ_,( + eZ + ,heH c + , (23) 
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where the Z± are the groups obtained by exponentiating the E a ,a G R±. As an example, 
the gaussian factor can be written as 

C_ = exp( zaE o)- (24) 

aER- 

When the representation matrix D(g) is applied to a greatest weight, the factor ( + is regarded 
as the identity. The z a can be taken as complex coordinates for the coset space G c /B + . It 
follows that the set of physically distinct states can have an one-to-one correspondence to 
the G c /B + . For a complex manifold we can choose complex coordinates on the manifold 
with holomorphic functions. And constructed on non-greatest weight states, coherent states 
are described by nonholomorphic functions of z. In the next section, after introducing 



supersymmetry on the compl 
mechanics in graphene 0, [l6] 



ex manifold in brief, we will describe supersymmetric quantum 



IV. SUPERSYMMETRY IN GRAPHENE 

Let us use the methods of Witten to introduce fermionic creation operators ip x ^ and ifj y ' 
which correspond to the differential forms dx and dy 0, 

dx <-> ip xj \0 >, dy <-> V^IO > . (25) 

In more detail ip x ^ performs the operation of exterior multiplication by da; while the adjoint, 
ip x , does that of interior multiplication by the vector dual to dx, say, d x 

i : d x -> ip x , S : dx -> ip x \ i : d y — > ip y , 8 : dy -> ip v ^. (26) 

The fermionic operators satisfy the anticommutation relations 

{^,^ f } = ^, tt,v = x,y. (27) 

On the basis of these definitions let us set up 

-(V>*-#s/) ="0*= 2^*' \{^x + i^ y ) = ^z = ]^P z - (28) 
while holding the Hermitian conjugate relations. 
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Now let us describe a supersymmetry over the complex Kahler manifold. On the complex 
manifold we take two supercharges 

Q 1 = d = ^ Z 8 Z , Q\ = 5 = -Vd- Z) 

Q 2 = d = ^d- Z , Q\ = 5 = -i) z d z . (29) 

Making use of anticommutation relations for fermions, these supercharge operators allows 
us to express 

85 + 58 = Q 1 Q\ + QlQi = -iv 2 , 

d5 + 5d = Q 2 Q\ + Q\Q 2 = -\v 2 . (30) 

Here it is easy to check that the cross terms such as 

85 + 58 = Q x Q\ + Q\Q x = 0, (31) 

do not have any contributions. 

In terms of the two supercharge operators, the Dirac operators D± in graphene is given 

by 

D ± = Tthv F {Q 1 + Q 2 ), = Q ± . (32) 

The Dirac operators can be described by means of the sum of the two supercharges which 
is the ordinary exterior derivative, d. And Q± become the supercharge of N = 1 super- 
symmetric quantum mechanics. It is not hard to check up that D± = Q± + Q± is given 
by 

t / d z \ 

Q± + Ql = T2ihv F (33) 
\d- z J 

which are equivalent to Eq. (6). Furthermore the square of Eq. (33) can lead to 

(Q± + Ql? = -4(^ F ) 2 f dz9 ~ Z ° ] = -h 2 v 2 F V 2 , (34) 

\ d- z d z ) 

where we have exploited Eq. (7). And hence in the sense of SUSY QM the Hamiltonian of 
graphene may be recapped in terms of the following form 

H ^ 2(Q ± + Ql) 2 = -2h 2 v 2 F V 2 , = ~^, (35) 
11 



provided that the mass, m* were defined by m* = in the last expression of the eq. (17). 
The Witten index is given by 

Index(d) = Tr ((-l) F e -' H ) (36) 

which accounts for the Euler number of the manifold as the exterior calculus of the de-Rham 
complex. In order to build up the supersymmety over the Dolbeault complex 22J, we need 
one of the supercharges, Q2 = 8. On the manifold of real dimension 2n, the index of the 
Dolbeault complex is given by 

Index(<9) = Tr f(-lf e -^ Q2+Q ^ 2 ) (37) 

This index is more interested in the the SUSY QM over twisted Dolbeault complex which 
is associated with deformation of the topology of the lattice on a graphene system . And 
in the next section we build the Atiyah-Singer and G index theorem as well as the energy 
eigenvalues on the deformation of the compact manifold. 



V. INDEX THEOREM IN GRAPHENE 



Let us describe the index theorem which gives an insight on the spectrum structure of 
certain operators such as the Dirac operators. In graphene this theorem enables us to have 
physical properties associated with the topology and geometry of the space in which the 
Dirac operators are defined. It provides the relationship between the analytic properties of 
the operator and the topological characteristics of the manifold. 



A. The Atiyah-Singer index theorem in graphene 

We illuminate the Atiyah-Singer index theorem by the method employed to the heat 
kernel expansion. The theorem furnishes a relation between zero eigenvalues of the Dirac 
operator of graphene and the total flux which goes through its surface. If the latter is 
connected to the genus of the surface through the Euler characteristic, we can find a close 
relation between the zero modes and the topology of the surface on a graphene system. 

Let us start with a Dirac operator given by 

KH° D, |. (») 
DO 

12 



Here D means an operator that maps a space M + onto a space M_ while is a map from 
M_ to M+. If D is an n x m matrix, becomes a m x n matrix. M + and M_ are the space 
of n and m dimensional vectors, respectively. Because we are focusing on the zero modes 
of K such that the solutions of the equations K\P = 0, let us define the number of different 
eigenstates of D with zero eigenvalue as r] + and the ones of as r\_. As a bookkeeping of 
chirality, the chirality operator 75 is defined as 



75 = • (39) 




Its eigenstates can have eigenvalue ± provided that they act on M±. 

In order to calculate the number of zero eigenstates in which we are interested, we take 
into account the operator K 2 that has the same number of zero modes as K. The K 2 can be 
given by a diagonal form 



K 2 = . (40) 




It is claimed that the operators D^D and DD^ get the same non-zero eigenvalues. To prove 
this statement, assume that DD^ijj = Xtp for eigenvalue A ^ 0. Then it follows that 

L>L>V = \ip — > D ] D{D ] i)) = A(L>V). (41) 

This means that the operator D^D gives the same eigenvalue, A, which corresponds to the 
eigenstate Dty. But it is not necessary to hold the case for A = when Dty might be zero 
by itself. 

Let us compute the trace of ^e"^ as followings: 

Tr( 75 e- tK2 ) = Tr(e~ tDtD ) - Tr(e~' DDt ) = 5>~ tA + - $>-* A -, (42) 

A+ A_ 

where A + and A_ indicate the eigenvalues of the operators D^D and DD\ respectively, 
and t is an arbitrary parameter. In the first step of the above procedures, 75 acts on the 
exponential so that it provides a +1 to the eigenvectors of D^D when they are placed in 
M + , and a —1 to the ones of DD^ when they belong to M_. In the last step the trace is 
evaluated by a sum over all the eigenvalues of the corresponding operators. Every non-zero 
eigenvalue of D^D is a one-to-one correspondence to an eigenvalue of DD^. Therefore all 
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paired terms of non-zero eigenvalues cancel out each other. There are left over the zero 
eigenvalues of each operators, resulting in 

Tr( 75 e- tK2 )=r ?+ -r / _. (43) 

In general we cannot determine difference between the number of zero modes. It is seen 
that the above result is independent of t owing to the cancelation of the non-zero eigenvalue 
term. 

Actually we should evaluate Index(K). In order to calculate it practically, we take an 
alternative method of heat expansion for calculating Tr(75e~* K2 ). It says that for general T 
and D on a two dimensional compact manifold we can expand 

Tr(te- t£) ) = ^- t J2th l (t,D), (44) 

l>0 

where Tr indicates the trace of matrices and the integration over coordinates of space, bi 
denote expansion coefficients. For f = 75 and D = K 2 , we have to return to an expression 
that is t-independent. For this t-independence, the expansion coefficients should vanish for 
all I except for I = 2 under the condition that all the t contributions is canceled out each 
other. This allows us to determine the coefficient b 2 from the first order term in t in the series 
of expansion. For the evaluation, let us take D as D = —ie^a u {y^ — ieA^). Here indicates 
the zweibein of curved surface metric g^ u that defines a local flat frame r] a p = e^e^g^ while 
is the Pauli matrix. denotes a gauge field. It follows that 

K 2 = -g^V,V u + \[Y,Y\F^ - \R (45) 

where R is curvature, and F^ v = d^A u — d v A^ means the field strength. V M indicates a 
covariant derivative with respect to gauge and reparametrization transformation. 
It is easy to see that the non-zero expansion coefficient b 2 is written by 

b 2 = Tr[ 7 5(^[7 M ,7l^ - \R)} = 2 I y B • dS, (46) 

where B is the magnetic field given by Bi = -e ^F^. The integration has been over the 
whole surface. These two independent ways for calculating Tr(7 5 e~* K2 ) allows u to arrive at 
the final formula of the index theorem 8j 

Index(K) = u + - p_ = J J B ■ dS. (47) 
14 



It states that the total flux which goes out of the surface is related to the number of zero 
modes of the K operator. The curvature doesn't contribute to the index formula because 
7 5 is a traceless operator. Particularly no contribution of curvature is to show an intrinsic 
property of two dimensional surfaces. The index theorem leads to an integer number on 
the compact surfaces. Therefore J J B ■ dS produces the total magnetic monopole charge 
in discrete values inside the surface under the Dirac quantization condition of magnetic 
monopoles. 

Using the index theorem, we observe the topological characteristics of graphene. Let us 
apply it to graphene and its geometric variants. In a certain configuration of graphene we 
have to account for the effective magnetic field on the graphene surface. Some plaquette 
deformations give rise to a specific circulation of the vector potential around a loop trajectory. 
Stokes's theorem helps us to have a relation between the circulation of the gauge potential 
around a loop 7$ and the flux of the corresponding magnetic field 

I A-df= f I B-dS, (48) 

where Si indicates the area. Hence going through the surface of graphene the total flux 
can be obtained from the fact that we can know the total number of deformations. It is 
necessary for us to have some information about the total number of plaguette deformations. 
The information is related to the topological properties of the surface through the Euler 
characteristics. Let us consider the Euler theorem in the following subsection. 

B. Euler theorem 

In general the Euler theorem gives rise to a relation between the structural information of 
a polyhedral lattice and its topological properties. There are a lot of proofs on this theorem. 
The most common methods are proofs on the basis of a reductions from the polyhedral 
lattice to the simpler one without changing its topological properties. 

Let us consider a lattice placing on a compact surface with a certain genus g. Then we 
can compute the number of deformations in a lattice necessary to create such a surface by 
applying the Euler characteristic. Let V, E and F be the number of vertices, edges and 
faces of the lattice, respectively, and N end , open ends. Then the Euler characteristic, x is 
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expressed by [8] 



x = V ^E + F = 2(l-g)-N end . 



(49) 



The second step of Eq. (49) is satisfied by the Euler theorem. It is easily to check that a 
single cut in the surface can have a reduction of the genus by one and increase the number 
of open ends by two, say, (g,N end ) —>-(<? — l,N end + 2), leaving the Euler characteristic x 
preserved. 

Let us apply the Euler theorem to the case of graphene molecules. There are three 
links on each vertex of graphene. Suppose that topological deformations such as pentagons 
or heptagons are present. Let us indicate the total number of pentagons, hexagons and 
heptagons by 71,5, and rij, respectively in the molecule. Then the total number of vertices 
is written by V = (5n^ + 6n 6 + 7n 7 )/3 when each fc-gon has k vertices and each vertex takes 
three polygons. Similarly, the total number of edges is expressed by E = (5n 5 + 6n 6 + 7n 7 )/2 
provided that each edge has two polygons. The total number of faces is equal to the sum of 
different polygons, F = n 5 + tiq + n-j. Combining these into the Euler characteristic, we see 
that 

n 5 -n 7 = 6 X = 12(1 - g) - 6N end . (50) 

This result reflects many facts. When equal numbers of pentagons and heptagons are in- 
serted, they do not make any change about the topology of the surface in the case that they 
cancel out. On a flat graphene sheet we can put two pentagons and two heptagons on it with- 
out changing the curvature of the molecule away from these deformations. This is consistent 
with the effective gauge flux approaches where pentagons and heptagons give opposite flux 
contributions. On the other hand, it is known that nontrivial topologies necessarily provide 
an imbalance between pentagons and heptagons. The genus zero configurations result in an 
excess of pentagons while high genus surface has an excess of heptagons. Genus one surfaces 
do not have any pentagons or heptagons at all provided that they are equivalent to a flat 
sheet. 

It is seen that Eq. (50) recaps the known result of a sphere with g = which leads to 
X = 2. This corresponds to a fact about 725 = 12 and nj = for the Ceo fullerene. For a torus 
with g = 1, we can have x = related to = n-j = in the case of the nanotubes. Thus 
no pentagons or heptagons may be required. If we account for the genus g = 2 surfaces, 
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then we can obtain x — ~ 2 where = and n-j = 12. In this situation equal numbers 
of pentagons and heptagons can be inserted without making any change of topology on the 
surface. 

Suppose that K is a Dirac Hamiltonian H. Then we compute the Index(H). The Euler 
characteristic term allows us to compute the gauge field term in the Index(H). It can be 
given by including additionally the contributions from the surplus of pentagons or heptagons. 
Hence the total flux of the effective gauge field can lead to 



, , „ . ... , r A • df 

2tt J J s 2vr ^ J 

= ~(n B - n 7 ) = 3(1 -g)- ^N end . (51) 

Z7T Z I 

The total number of zero modes is equivalent to the sum contributed from each subsector of 
a Dirac operator. As a consequence, by adding the two contributions, we arrive at the index 
of the Dirac Hamiltonian that describes the graphene molecule. The Index(H) is expressed 
by 

Index(H) = v+ - v_ = 6(1 - g) - 3N end (52) 

which is consistent with the exact number of the zero modes if v + = or z/_ = 0. Therefore 
we have obtained the theorem which relates the number of zero modes existed in a certain 
graphene molecule to the topological characteristics of its surface. 

This result provides the number of zero modes for the familiar cases of graphene molecules. 
For example, since a fullerene takes g = and N en d = 0, it is expected that it has six zero 
modes which correspond to the two triplets of Cqq and of similar large molecules. For the 
case of nanotubes, we have g = and N en d = 0. This results in y+ — z/_ = which is in 



agreement with previous theoretical and experimental results [20|, |2lj]. The index theorem 
gives rise to a surprising relation between the topology and the presence of magnetic flux 
which is effectively inserted in graphene molecules by geometrical deformations. The number 
of these deformations can be associated with the general topological characteristics of the 
lattice surface. They are related to the zero modes of a general graphene molecule with the 
genus and the number of open faces of its surface. 
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C. G index theorem 



To build up supersymmetry over the Dolbeault complex which are associated with de- 



formation on topology of the lattice 22j, we need one of the supercharges, Qi = 8. The 
index theorem leads to a topological invariant under deformations on Q 2 and plays an es- 
sential role in formulating the SUSY QM over the twisted Dolbeault complex caused by the 
deformation in a graphene system. 

Let us describe the general statement of the G-index theorem. Suppose that Q and 
are supercharges which have a map from a space of bosonic states to a space of fermionic 
states and vice versa. Further let us take a Lie group G generated by Gi. Gi satisfy the 



commutation relations with Q and Q} such that [17] 



[d, Gj] = z/*G fc , [Q, G { ] = 0, [Qt, d] = 0, (53) 

/j* are structure constants. Gi has also a commutation relation with the Hamiltonian of Eq. 

(5) 

[Gi,H] = 0. (54) 
A supercharacter can be expressed as a supertrace of the group elements 

E(e i6iGi ) = Tr ( (-if e^V 11 ) . (55) 



The ordinary trace requires us to take a limit of t — > 0. But Eq. (55) does not have 
any dependence on t because the non-zero energy levels can not make contributions due to 
canceling in pairs between the bosonic and fermionic sectors. But the zero-energy levels can 
only give contributions. And these are not dependent on t. 

The character is topologically invariant under deformations of the operators Q and QL 
On deformation let us consider the supersymmetry generator Q of the ordinary N = 1 SUSY 
QM on a manifold. In locally geodesic coordinates Q is expressed by 

Q = d = V^fy. (56) 

When deforming Q, it is changed into the new operator 

Q s = d + si K = V ] d» + sJC^, (57) 
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where K, denotes a Killing vector field. And then the deformed Hamiltonian is given by 

PL = 2(Q S + Q\) 2 , 

= ~d 2 + s 2 \IC\ 2 - l -s[^l]d^K, v - l -s^^ u ]d v K,. (58) 

Now suppose that we decompose the deformed N = 1 SUSY operator into holomorphic 
and antiholomorphic sectors because of d = d + d. Then the G-index theorem can be 
associated with a high degree of symmetry on graphene molecule. Let us take a Lie group 
G generated by Gi which satisfies the commutation relations with Q 2 and Q\ such that 

[G h Gj] = i/*G ft , [Q 2 , Gi] = 0, [Ql d} = (59) 

where mean structure constants. On graphene, there exists the Hamiltonian H such that 

H ee -2h 2 vl(Q 2 + Ql) 2 , [G h H] = 0. (60) 

On deforming Q%, this can be changed into a new operator 

0,2s = ^~ z d- z + s/C%-, (61) 

where s is a real parameter. K, denotes a Killing vector field associated with the deformation. 
In the complex coordinates z and z, the Killing vector field is given by K = —yd x + xd y = 
i(zd z - zdg). 

Under the deformation let us consider a graphene system corresponding to the K + . Then 
we can obtain a deformed Hamiltonian H +)S on the Dolbeault complex. For convenience, 
after replacing s by is, the deformed Hamiltonian is given by 

H +iis = -2h 2 vUQ+,2is + Q{,2 is ) 2 

= -2h 2 v 2 F (-2d 2 Zz + \s\ 2 \z\ 2 ) - 2h 2 v\ [s\^~ z , + s(zd s - zd z )) . (62) 

This operator gives rise to eigenvalues j^l, Q 



E+ : nml — —2H 2 V 2 ? 



i r 

n+-) + (m+- ) +.s-/ + .s(±i; 



(63) 



Depending on the choice of s > or s < 0, we can consider two cases for zero eigenvalues 
satisfying the topological invariance imposed by the index theorem. First if s > 0, we should 
take —1 for the [ip* x ,tl)g]. Now we get in the bosonic sector. We can have zero eigenvalues 
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in the case of I = 0, —1, —2, • • • . Second if s < 0, we have to choose +1 for [ip^ z , 4>z\- And 
then we can get zero eigenvalues for I — 0, 1, 2, • • • . 

According to the choice of s, let us take two Hamiltonians is = 2 isQ+,2is an d 
H+ is = Q+,2is<5+ 2is as t wo superpartners which imply the Z 2 grading over the Hilbert space. 
Here we have expressed Q + ^is = —iV2frvF(Q+,2is + Q+2is)- Then we want to calculate the 
eigenvalues for the up-spin and down-spin components 

K\is\^ + \nml >= E\\ nml \^ + \ nml >, (64) 

where E^ n+lml > E^ nml > E^ 0Q0 . Suppose that E\ Q00 is zero. Then making use of the 
relations 



Q+,i2s^+,is — Q+^isQ+^isQ+^is — H+ t i s Q+,2is, (65) 



we can have 



H +,i,Q+,2w|V'+,„ fn l > = §+,2wHj_ ) JV>£ )nfnl > 
= El,n m lQ+,2isWl,nml > • (66) 

This means that if nml ^ 0, Q+ : 2is\4>\. nm i > 1S an eigenstate of H^j_ is . And similarly 

^+,isQ +,2is\' l f ) +,nml > = Q+,2is^+,is \^+,nml > 

= Ei jnml Ql >2is \^ +inml > . (67) 

If E^ nml 7^ 0, Q+ i 2i S |V ; + nmi > * s an eigenstate of H^ is . Therefore, for non-zero eigenvalues 
there exist the up-spin and down-spin eigenstates in pair. They form a supermultiplet 
connected by the supercharge Q+$is- The up-spin (down-spin) sector may be described as 
the bosonic (fermionic) sector at the K + point. 

In the case of the zero eigenvalues for the down-spin sector, we should investigate both 
E+ ooo = an d -^+000 7^ separately. Let us assume that E\ 000 ^ 0. Then the lowest 
eigenstate of H^ s is expressed as |^+ j0 oo > ^ < 5+,2«s|V'+ i m i > • Hence, there exists a su- 
permultiplet between the states l^+n+im; > an d IViw > ' which have the same energy 
E^ n+lml = E^ nml , and we obtain energy eigenstates |9| 

\^+,nml > = I = =Q+,2is\lpX tn+lrn i >, 

e\ 1 



+,n+lm£ 

',t 



\r + ,n + lml > = I , Q+,2is\n,nml > ( 68 ) 

\e\ 



,nml I 
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for n > 0. If 000 = and Q' + 2 isW+ ooo >= 0, the relationships are written by 

1 ^ 

Y \ E +,nml\ 

\Anml > = A Qi^inml > (69) 



,nml I 

for n > 1. Similarly, one can repeat eigenvalue problem for H-^, which corresponds to 
the X_ point. As a relation between the up-spin and down-spin eigenstates, the bosonic 
(fermionic) sector is regarded as the down-spin (up-spin) sector at the K_ point. And hence 
there exists the 4-fold degenerate energy spectrum. 

VI. DEFORMED ENERGY EIGENVALUES 

AND UNCONVENTIONAL QUANTUM HALL EFFECT 



On a sheet of graphene, let us consider the problem of magnetic field concentrated on 

ch_ 
eB' 



a thin cylindrical shell of small, but finite radius Ib = yf|- The corresponding vector 
potential is given by a = (—y,x)/2l B on the two dimensional plane of graphene. The 
problem in question is to compute the eigenvalues of the Dirac Hamiltonian in the field of a 



fractional magnetic flux on the graphene sheet 23|, |24|. Now under the fractional magnetic 



flux, the eigenvalues for n = 0, 1 = 0, m = are given by 



/4,qoo = V<ooo = 0, (70) 

and 



^l,n+100 = V ^l.nOO = V E -, n+100 = V E ~, nOO = ±hw iB Vn + 1 (71) 

for n > 0, 1 = 0, m = 0. Here W\ B = V ^' F . Equation (70) tells us that there is one zero-energy 
state only in the case of up-spin fermions but not in the case of down-spin fermions at the 
K + point. At K_ point we can have one zero- mode state for down-spin fermions but not 
for up-spin fermions. The magnetic field direction at K + is opposite to that at K_. The 



,8 



zero-energy state may have the four-fold degeneracy emerging from electrons and holes 
Since the LL of the zero-energy states becomes half-filled, no one would observe plateau 
at v = 0. But by index theorem, the flux quanta produce 4r (r = 0, 1, 2, • • • ) zero-energy 
states. The 2r states of these are occupied. The flux quanta lift the 2r states to the Fermi 
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energy. And then they can be removed by doping. The degeneracy between electrons and 
holes would be removed. We could observe the Hall plateau at u — because holes are 
occupied before electrons. And hence we can describe an experimental observation of the 



Hall plateau emerging at v — 0. 



id 



On the basis of the index theorem, we compute the energy spectrum of the deformed 
Hamiltonian, H +jis = —2h 2 Vp(Q +j 2is + Q+2is) 2 a t the K + point. The up-spin states of zero 
energy are |0 >, |1 >, |2 >, • ■ ■ — 1 >. They are degenerate in \^\ 0m i >■ On the 
other hand, for down-spin states, we have to describe two cases. As the first case, assume 
that = 0. Then, there do not exist any zero-energy states. So we may construct the 
supermultiplet given by Eq. (68). In the other case, if ^ 0, the zero-energy states are 
given by |0 >, |1 >, |2 >, • • ■ , \j^ — 1 > as degenerate states of |^+ 0mi > • Therefore, these 
result in the (j^ + j^)— fold degeneracy in the zero-energy state for fermions at the K + point. 
This degeneracy implies the exact correspondence between p fermions and p fermions under 
deformation. Similarly, we can investigate the energy spectrum of the zero-energy states at 
the K_ point. 

In order to generate the up-spin and down-spin states of zero energy the deformed su- 
peroperators are written, in terms of the original supercharges, by 

Q+,2is = Q+,2isQ+,2isi Q\,2is = Q+,2isQ+,2isi (^2) 

where j 1 " and are integers such as j 1 " > jK In terms of the Q+ 2 is the state \n > is given 
by 

In >= -U(QU) n |° > • (73) 
vn! 

Now, let us calculate the deformed eigenvalues by using Eq. (73) and solving the eigen- 
value problems of Eq. (64) for the up-spin and down-spin components. In the bosonic sector, 
the deformed energy eigenvalues are expressed by 



'n.ooo = V<ooo = 0, (74) 
for n = 0, 1 = 0, m = 0. And we can have 



K.nml = \l E -,nml = ±fc«>lfl \ / [(7l-l)!]2 ^ > 
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for n > 1,1 — —(y — l),m — 2(p — 1). In the case of the fermionic sector, if y = 0, the 
deformed eigenvalues are 



K,nml = \/ E -,nml = ±hw l B \j ^yjp (76) 



for n > 0, / = y, m = 2y. If y 7^ 0, the eigenvalues are given by 



'4,000 = V^-,000 = °' ( 77 ) 
for n = 0, I = 0, m = 0. And we obtain 



(n + j± - l)!(n + jt - 1)! 

-,n+lml - [(„_ 1)!]2 



for n > 1,/ = — l,m = 2(j^ — 1). We can check up that there exists A(j* + j^-fold 
degeneracy in the zero-energy states and 4-fold one in all other states. 

Among the energy spectrum given above, let us account for the special cases of = 1, 
and y = and = 2, and y — 0. In the case of y* — 1, and y = it is not hard to check 
up that the energy eigenvalues are given by Eqs. (77) and (78). These results correspond 
to the spectrum of the monolayer graphene. For the case of y — 2, and = 0, the energy 
spectra are expressed in terms of 



£[,000 = \/<ooo = 0, (79) 



for n = 0, 1 = 0, m = 0. And we can obtain 18] 



El,nml = V E -,nml = ±^B >/n(n + 1) (80) 

for n > 1, / = —1, m = 2, while having 



*W = V^lw == ± hw iB V(n + l)(n + 2) (81) 



lene affected 



for n > 0, 1 = 0, m = 0. These energy spectra are eigenvalues of the bilayer gra 
by the deformation, and are in agreement with the results in the literature [10 
the basis of the index theorem we have shown that there exist the j'-fold and j^-degeneracy 
in the zero-energy state at the K + point and similarly at the K_ . And hence we can obtain 
the QHE characterized by [9J] 

a xy = v^,v = ±4(|n| + ^). (82) 
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VII. SUMMARY AND CONCLUSION 



We presented the electronic properties of massless Dirac fermions characterized by geom- 
etry and topology on a graphene sheet in this chapter. Topological effects can be elegantly 
described by the Atiyah-Singer index theorem. It provides a topological invariant under 
deformations on the Dirac operator and plays an essential role in formulating supersym- 
metric quantum mechanics over twisted Dolbeault complex associated with the deformation 
on the topology of the lattice in a graphene system. Exploiting the G-index theorem and 
a high degree of symmetry, we explained deformed energy eigenvalues in graphene. The 
Dirac fermions result in SU(4) symmetry emerging out of both the pseudospin and spin as a 
high degree of symmetry in the noninteracting Hamiltonian of monolayer graphene. Under 
the topological deformation the zero-energy states emerge naturally without the Zeeman 
splitting at the Fermi points in the graphene sheet. Thus we observed an emergence of 
a higher degree of hidden symmetry under the topological deformation in graphene while 
the pseudospin is a good symmetry at the K and K' points in graphene. In the particular 
SU(2) of the pseudospin, the SU(2) is the exact spin symmetry of each Landau level. In the 
case of nonzero energy, the up-spin and down-spin states have the exact high symmetries 
of spin, forming the pseudospin singlet pairing. The pseudospin can play a key role on the 
physics of the n = LL in the graphene sheet. The valley pseudospin degeneracy can lift 
only at the zeroth LL. The 4-fold degeneracy can be removed in the zero-energy states of 
monolayer graphene. If the mass terms were taken into account, the four-fold degeneracy 
can be removed in the zero-energy state of monolayer graphene. We can exploit this to 
understand the emergence of a Hall plateau at n = in the experimental observations. But 
the four-fold degeneracy is not removed in the higher LLs. Including the Coulomb interac- 
tion, we can lift the degeneracy. The pseudospin symmetry SU(2) is broken to U(l) x Z 2 . 
Therefore the total symmetry gives rise to SU(2) sp j n x (U(l) x 7i) psuedospin while the spin 
symmetry SU(2) remains to be exact. Hence we understood the peculiar and unconven- 
tional quantum Hall effects of the n = Landau level in monolayer graphene on the basis 
of the index theorem and the high degree of symmetry under the topological deformation 
without the Zeeman splitting. It would be very interesting and quite possible to apply the 
present approach to investigation of the composite Dirac fermions and fractional quantum 
Hall effects in graphene. 
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